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Abstract 

We introduce and study a simple unbalanced holographic superconductor model with two 
scalar order parameters. The attention is focused on the possibility of coexisting orderings 
corresponding to the concomitant condensation of two scalar operators. Through a probe 
analysis we show that an attractive or repulsive direct interaction between the two bulk 
scalars leads respectively to competition and enhancement of the associated condensates. 
The system at hand is a toy model for studying generic multiple ordering in a strongly 
coupled context and some comments are given about its applicability to the ferromagnetic 
unconventional superconductors. 
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1 Introduction 

The AdS /CFT correspondence offers us the possibility of studying quantum field theory 
holographically, that is, in a dual perspective involving an appropriate gravity theory on 
a spacetime with higher dimensionality. According to the AdS /CFT conjecture, a classi- 
cal gravitational theory on an asymptotically AdS vacuum can be employed to compute 
correlation functions of the dual quantum field theory thought of as "living" on the con- 
formal boundary of the asymptotic AdS* geometry [lj|2]. This is true in an appropriate 
dynamical regime involving small curvature and weak coupling on the gravity side and 
strong coupling in the dual quantum field theory. 

Holography, as a framework to access the strongly coupled regime of quantum field 
theory, has received vast attention since it provides us with a new approach to study 
interesting toy models for real physical systems whose dynamics happens to be strongly 
coupled. These strongly coupled models are significant both in condensed matter and in 
QCD. 

A particularly important field of application for holographic analysis is given by the 
physics associated with quantum phase transitions. Within this context, models for high- 
T c superconductors have been widely studied holographically. Their non-standard prop- 
erties (i.e. not describable within the BCS framework) are thought of as being possibly 
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intimately related to the collective and strongly coupled dynamics of the system. The 
first model for a holographic superconductor has been introduced in [3j[4], see [5j[6] for 
a wider perspective. The ultimate and ambitious aim of the holographic approach to 
unconventional superconductors would be to shed light on the high-T c mechanism of su- 
perconduction or at least provide us with useful effective models. Indeed, already on a 
phenomenologically oriented level, holographic studies offer the possibility of handling 
quantitatively some macroscopic properties of strongly coupled superconducting systems, 
such as the thermodynamic and transport properties. 

As studied in it is possible to extend the minimal holographic superconductor 
[4|[8] and introduce an additional gauge field in the gravitational theory. According to 
the holographic dictionary, such new field represents an effective way to account for a 
second chemical potential in the boundary theory. Indeed, the extended holographic 
model introduced in u\ is meant to study the strongly coupled physics of unbalanced 
mixtures where more than one (in that case two) chemical species are related to different 
chemical potentials^ 

Unbalanced mixtures are important both in the condensed matter context and in 
QCD. In the latter the two bulk Abelian gauge fields correspond to two boundary global 
£7(1) symmetries which can describe the baryon number and isospin symmetries. Albeit 
in this paper we are going to focus the attention (and the language) on the condensed 
matter context, it should be kept in the back of one's mind that the holographic toy model 
under study allows for a wider range of applications. 

According to the unbalanced superconductor interpretation of [7], the two chemical 
potentials are read as the sum and the difference of the chemical potentials for spin-up 
and spin-down electrons. In this sense, the first bulk gauge field is associated to a global 
"electric" symmetry while the second bulk gauge field accounts for an effective "spin" 
symmetry of the boundary theory. A non- vanishing difference among the chemical poten- 
tials for spin-up and spin-down electrons is what makes the system unbalanced. As pointed 



out in 11 , in electronic systems at low energy the spin-orbit coupling is suppressed lead- 
ing to an effective decoupling between spin and spacetime rotations. Therefore, even in 
a context where spatial rotations happen to be broken, the low-energy SU(2) spin sym- 
metry is a valid approximate symmetry of the IR system. The "magnetic" U(l) in our 
holographic model can be thought of as, for instance, the 03 component of the SU(2) spin 
symmetrjj^J 

Unbalanced systems are particularly interesting in view of studying mixed spin-electric 
(in one word "spintronic" ) transport properties at strong coupling and the possible emer- 
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1 A seminal paper employing a holographic approach to the study of unbalanced systems is 
2 In the present treatment the ferromagnetic-like order parameter is effectively described with a scalar, 

namely, without loss of generality, we understand the choice of a fixed spatial direction with respect to 

which spins arc projected. 
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gence of inhomogeneous FFLO phases where the superconducting condensate acquires 
a non-trivial spatial modulation. FFLO phases at weak coupling were introduced and 
theoretically studied by Fulde, Ferrel, Larkin and Ovchinnikov [18] . The first suggestion 



about the possibility of FFLO-like phases in holography has been advanced in 25 ; further 
holographic studies have been carried out in [7] and separately in 26 , 27 . 



1.1 Motivations for a second order parameter 

The purpose of the present paper is to generalize the minimal setup for an unbalanced 
holographic superconductor adding a second scalar field to the model introduced in [7]. 
A scalar operator which undergoes condensation describes the occurrence of a non-trivial 
order parameter. For instance in the minimal holographic setup the charged scalar breaks, 
upon condensing, the electric U(l) symmetry and it is naturally interpreted as a super- 
conducting order parameteij^J In the presence of a second scalar we access the possibility 
of a more complicated phase diagram where two orderings can occur. The dynamics and 
mutual relation between the two condensates is one of the most interesting property to be 
studied; in other words, one would like to investigate whether the two orderings compete 
or enhance one another. 

The point is so far very general. Before making it more specific it is important to 
underline the versatility of a toy model possessing two order parameters. By simple 
changes of the parameters (e.g. the charges and masses of the two scalars) one is able to 
describe different and interesting physical situations. Indeed, the toy model at hand could 
admit various physical interpretations and allows us to address the problem of ordering 
coexistence at strong coupling in a rather general fashion. 

The focus of the present treatment is nevertheless on a specific case particularly rel- 
evant for condensed matter applications: We take the first scalar ip to be charged under 
the "electric" U(l) and neutral to the "magnetic" U(l) while for the second scalar A we 
consider the opposite charge assignment. Hence ip is meant to account for an electrically 
charged and magnetically neutral condensate, like for instance a condensate arising from 
an s-wave Cooper-like pairing. On the other hand, A could represent electrically neutral 
magnetization)^} 

A pivotal question in the experimental and theoretical research on high-T c supercon- 
ductors has been the interplay between superconductivity and magnetic orderings (see 



for instance 22 and references therein). From the theoretical standpoint, standard ef- 



This superconductivity claim can be supported with direct analysis of the DC diverging conductivity 



below T c , see for instance 13 for further details. 

4 Notice that we have just introduced a "magnetically charged" order parameter accounting for 
ferromagnetic-like ordering. One could also consider a neutral order parameter with the aim of describing 
antiferromagnetic-like orderings, see [9l[TTl . 
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fective techniques d la Landau- Ginzburg have usually been employed as, for instance, in 
the perturbative study of spin-density-waves in superconducting cuprates 2~Ip Comple- 



mentarily, we adopt here an effective and minimal approach which nevertheless describes 
an intrinsically strongly coupled medium which allows for coexistence of multiple order 
parameters. 

The general standard picture arising at weak coupling tends to discourage the concomi- 
tant occurrence of superconducting and long-range ferromagnetic order [24] . Moreover, in 
the weak coupling analysis, standard Fermi liquids present usually few instabilities and, 
more importantly here, the occurrence of an instability tends to disfavor further instabil- 
ity. This is intuitively due to the fact that when the Fermi liquid undergoes an instability, 
a mass gap is usually produced. Such mass gap tends to prevent, or at least discourage, 
further instabilities (9l. 

In a strongly correlated system however the picture might change radically. In partic- 
ular, in the panorama of high-T c superconductors a wealth of coexisting orderings such 
as magnetic orderings or striped phases are possible. The recent discovery of materials 



where itinerant ferromagnetism and superconductivity could be even cooperative (see 14 
and references therein for details) is particularly interesting. 



2 The minimal holographic unbalanced superconduc- 
tor 

We start introducing the holographic model for the minimal unbalanced superconductor 
in 2 + 1 dimensions in a rather succinct way; we refer to |7| for any further detail. The 
gravitational 3 + 1 dimensional bulk action is: 



S 



2k\ 



dx \f—g 



6 

L 2 



1 

4' 



R ~t~ ~7t -F nh F ab - -Y nh Y° 



— m% } ij>*ip — \dip — iqAAipy 



(2.1) 



where we have an ordinary Abelian Higgs model with gauge field A (F = dA) in the 
presence of a negative cosmological constant —6/L 2 plus a second gauge field B whose 
field strength is Y = dB. We underline that the scalar ip is minimally coupled to A and 
uncharged with respect to B. 

We study the model considering the standard radial ansatz: 



ds 2 = -g(r)e- xir) dt 2 + —(dx 2 + dy 2 ) + 

L 9{r) 



(2.2) 



5 A study of a Landau-Ginzburg approach in the case of two order parameters is analyzed in [10] . In 
addition, the paper contains a list of interesting physical applications of the general problem of coexistence 
of multiple orderings. 
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ip = ifj(r) , A a dx a = <f)(r)dt , B a dx a = v(r)dt , (2.3) 



and solve the equations of motion descending from (2.1). We skip all the analysis which 
has already been described in [7] and focus on some interesting results regarding the 
relation between the imbalance and the condensation of i/j. 

We remind the reader that the boundary values for the gauge fields A and B are 
respectively interpreted as the overall chemical potential fi and the chemical potential 
imbalance 5fi of the two fermion families (namely "spin- up" and "spin-down" electrons). 
The thermodynamics of the finite temperature boundary field theory is mapped into 
the thermodynamics of the dual black hole and vice versa [29]. To simplify the following 
expressions we henceforth consider the AdSs + i curvature radius to be equal to one, namely 
L — 1; this can be done in full generality and amounts to a choice of length unit. 



2.1 Phase diagram: the imbalance hinders the scalar condensa- 
tion 

The phase diagram of the unbalanced holographic superconductor in the case of = — 2 
is plotted in Figure [TJ As noted in [7] (where the diagram was obtained) it is possible to 
observe that, for this specific mass value, scalar condensation appears for any value of the 
imbalance. In other terms there is no Chandrasekhar-Clogston bound [19] or, equivalently, 
for any Sfi/fj, there is a finite critical temperature below which the system undergoes the 
superconducting transition. Note however that the bigger the imbalance, the lower the 
critical temperature and therefore a larger value for Sfi/fi discourages condensation. This 
observation is in line with the interpretation of the condensate as formed by Cooper- 
like spin-singlet pairs and indeed it matches with the weak coupling expectation (see for 
instance |20~]). 

The instability leading to scalar hair formation is associated to the violation of the 
Breitenlohner-Freedman bound for a scalar field living on the near-horizon AdSi geometry 
[l~5] . The study of the near horizon geometry of the doubly charged Reissner-Nordstrom 
black hole solution at T = leads to the following condition for instability}^] 

rhl < ~, (2.4) 
where m 2 is the effective IR mass of the scalar field given by 



™>l = m % ~ 7. 2q la, ' ( 2 - 5 ) 



3 Some details are given in Subsection 4.1 of [7]. 
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Figure 1: On the left, phase diagram of the unbalanced superconductor; the shaded region 
corresponds to the superconducting condensed phase separated from the normal region 
by a line of second order phase transitions. Temperature is normalized with respect to 
the balanced case, i.e. T c (5fi = 0). On the right we have a plot showing the difference 
between the square of the IR effective mass and the bare mass of the bulk scalar field; 
the plot corresponds to q = 2 and fj, — 1. 



plotted in Figure pi Since the model we consider has 



m*p = —2, inequality (2.4) is always 
satisfied and, as already noted, the RN doubly charged black hole is always unstable 
towards hair formation at sufficiently low T. Nevertheless, let us underline that for 
larger values of the chemical imbalance the effective mass (2.5) is larger as well; then, 
increasing 5fi/fi the condition for instability, though always satisfied, is met with a smaller 
margin. This intuitively leads one to think that the T = doubly charged RN is "less 
unstable" when the imbalance increases and accordingly the condensation requires a lower 
temperature. 



The holographic unbalanced superconductor with 
two order parameters 



We introduce a second scalar field A in the gravitational model (2.1). It is minimally 
coupled to the gauge field B and neutral with respect to A; in addition, we consider a 
direct coupling 7 between the two scalars. In practice we add the following terms to the 



action (2.1 ) 



S x = — I dx A yf^ [-mlX^X - \dX - iq B B\\ 2 - 7 ^ip A f A] . (3.1) 



We take m\ = —2 so, in the specific case at hand, if} and A have the same mass value; 
this value is in agreement with many instances where one studies the AdS^/ 'CFT 3 cor- 
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respondence in a UV completed framework, namely considering string theory consistent 
truncations 30,31 . In accordance with the ansatz (2.3) we choose a radial ansatz for A 



as well, 

A = A(r) . (3.2) 

Putting 7 = we simply recover a "double-copy" of the standard superconductor jij[8] . 
In other terms, we have two scalar/gauge vector sectors, namely (if), A) and (A, B), both 
minimally coupled with gravity. Let us observe that gravity couples the two sectors but, 
as far as this paper is concerned, we will neglect the gravitational interaction (i.e. we 



work in the probe approximation, see Section 3.1) and the direct coupling 7 turns out to 
be cruciaQ 

According to the "magnetic" interpretation of U(1)b, the A condensate describes a 
magnetic order parameter. Hence, the coexistence of the two condensates, if) and A, 
corresponds to a ferromagnetic superconducting phase^J Without pretending to be par- 
ticularly close to any specific realistic model, let us however observe that the "magnetic" 
degrees of freedom giving rise to a condensate resemble itinerant ferromagnetism. Said 
otherwise, the degrees of freedom that we describe in terms of the "magnetic" condensate 
are not localized or fixed to any spatial pattern or lattice; in this respect, recall that our 
ansatz does not present any specific spatial features, neither periodicity nor lattice-like 
structures^} A condensate is likely to describe a plasma of degrees of freedom rather than 
an array of localized spins; this can be thought of in analogy with the fact that if) is 
interpreted as a condensate of (clearly itinerant) Cooper-like pairs at strong coupling. 



The total action is obtained summing (2.1) and (3.1) 



Stot — S + S\ . (3-3) 
The equations of motion descending from S tot are 

0"+(| + ^0 / -^^ 2 = O (3.4) 



r + (ff:_£ + HV + ^ _ !!i±^! v , = o ,3.5) 

9 2 rj g 2 g 



7 It would be interesting two consider instead, or in addition, also a direct coupling between the two 
gauge fields. We defer this to future work. 

8 A treatment of ferromagnetic superconductors performed by means of Hubbard-like Hamiltonians is 
described in [14] ; we refer to this paper and references therein also for a list of interesting and recently 
investigated ferromagnetic superconducting systems. 

9 Any phenomenon of commensurate magnetism appears then to be outside of the present treatment. 
Indeed, with "commensurate" we mean that the spatial modulations of the order parameter correspond 
to the periodicity of an underlying lattice. However, there is an interesting interplay between possi- 
ble modulated s-wave superconducting phases and incommensurate spin-density-wave ordering (see for 
instance [23]) which could be worth studying from a holographic viewpoint. 
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/ , ,/2 , x/2 A rq\4> 2 i) 2 e x rq B v 2 X 2 e x 
X + rip + rX -\ 5 1 t. = 



y + v- + ^ 



/2 e X y/2 eX 



+ 



+ - + - 2 

gr r z 



rn 



I /2 , 7 x2 ,2 , m L 2 , gi^ 2 2 e x . q 2 B X 2 v 2 e x 



+ — ^ + —X A ip + —A + 



+ 



2^ 2 



A" + 



fli 2 r 



g- 
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v" + 



X' 2 
— + - 
2 r 



v '_ 2 A X 2 v = 
9 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



To solve them we consider the ansatz (2.3) and (3.2). We study black hole solutions and 
the horizon is defined by g{rh) = 0. The equations of motion have the following scaling 
symmetry 



r -> ar , (t,x,y) -> -(t,x,y) 

a 



9 -» a 



— >■ a< 



(3.10) 



In an asymptotically AdS^ background, the near boundary behavior of the gauge and 
scalar fields is 



P>A 



PA 



+ ... , v = [i B 



Pb . , Ci C 2 « 

— + ..., ip = — + — + ... , A 



* + ^ + ... (8.113 



To underline the generality and symmetry of the present model, we have just adopted a 
less specific notation trading respectively p, p, 5p and Sp (used in the previous Sections) 
for p A , Pa, Pb and p B . 

Note that, under the scaling ( 3.10[ ), the near-boundary quantities behave as follows: 

Pa,b -> a/i^s , Pa,b -> °? Pa,b , (Ci, A) -> a(Ci, D x ) , (C 2 , £> 2 ) -> a 2 (Ci, -Di) • 

(3.12) 

We will be concerned in studying the spontaneous condensation of the scalars and we 
define the condensates as follows 



O, 



V x = V2 D x 



(3.13) 



where x = 1,2 and the conventional factor \/2 is introduced to comply with the existing 
literature; the dimension of the two kinds of condensates is \mass] x . 
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Figure 2: Examples of coexisting ip (black) and A (gray) condensates plotted by performing 
a sweep in temperature at constant charge densities. On the left/right plot the UV 
subleading/leading terms for the scalars have been respectively fixed to zero. The A 
(gray) plot is interrupted at a temperature value which is below its own condensation 
temperature because our sweeps are technically performed "moving" ^>(r#), i.e. the 
value of ij} at the horizon; indeed, we have ip(rjj) — at and above the ip condensation 



temperature (see Section 4.2). The temperature is here normalized with respect to the 



critical temperature for ip condensation. 



3.1 Probe analysis 

We do not have analytic solutions for the system of equations of motion, then we tackle 
the problem numerically. We work in the probe approximation that consists in neglecting 
the effect of the gauge and scalar fields on the geometry. Such approximation becomes 
exact in the limit of large charges q^, qs — > oo with qx^, qAA, g_eA, qsB kept fixed 
and finite. Indeed, in this limit, the terms involving the gauge and scalar fields inside the 
equations of motion for the metric become negligible pi; the metric is therefore a fixed 
background describing an uncharged AdS'-Schwarzschild black hole, namely 

g{r) = r 2 [l- r ^j , x(r) = 0, r H = %xT , (3.14) 

where the temperature T is linearly related to the black hole horizon radius r#. 

The equations of motion of the two gauge fields and the two scalars are ordinary second 
order differential equations, namely 

0" + %' _ _M ^ = v „ + 2 v , _ 2g| x2v = (3 15) 

r g r g 

r+ ( 2 - + l)^ + *l^-^^i- = (3.16) 
V r 9 J 9 9 



9 



X"+( 2 - + 9 -) \' + q jv*\- ml + ^\ = (3.17) 
\r 9 J 9 2 9 

Albeit the space of solutions of the system is spanned by eight parameters, we need 
to impose some constraints dictated by consistency reasons. According to the definition 
of black hole horizon we have g(rn) = and then, in order to avoid unphysical divergent 
quantities we must impose A(th) = B(th) — pj. Smoothness of solutions at the 
horizon implies a further consistency constraint [6): multiplying the scalar equations by 
g(r) and then considering them at r = r# one obtains that the value of a scalar field and 
its derivative at r H are not independent quantities. 

All in all, we have a four-parameter family of consistent solutions; as the consistency 
requirements impose some conditions at the horizon, our numerical code is build to propa- 
gate solutions from the horizon towards the boundary. The four "free" inputs are therefore 
the following horizon quantities, 

<f>'{r h ) , v'{r h ) , \{r h ) , $(r h ) . (3.18) 

Now we turn the attention to some physical features. In our analysis we always require 
that either the leading or the subleading term in the UV expansion of both bulk scalars 
is set to zero. We choose m 2 = —2 for both scalars, this implies that we fall inside the 
window where two quantizations are allowed [12] ; on practical grounds this means that 
we have the possibility of interpreting either the leading fall-off as the source and the 
subleading as the VEV of the corresponding dual operator or vice versa. Even without 
invoking any more sophisticated argument, requiring zero scalar sources corresponds to 
ask for unsourced operators VEV's, i.e. spontaneous condensations. Furthermore, we 
decide to work in the canonical ensemble and therefore the charge densities (i.e. the 



subleading terms in the UV behavior of the gauge fields, see Eq. (3.11 )) are imposed from 



outside. The four parameters (3.18) are therefore fixed and the code implements such 
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physical inputs by means of a shooting methoc 

An important observation is in order: As already noted, because of the scaling (3.10), 
the physically meaningful quantities need to be scaling invariant. When solving the 
holographic system we want to study a specified thermodynamical state (pa,Pb,T) but 



we are insensitive to an underlying scaling factor acting as in (3.10). Such feature can 
be read as emerging from the underlying conformal structure of the strongly coupled 
quantum field theory. 

We define the following "physical" (i.e. scaling invariant) quantities 

T O V 

T= —L O x = - - Ux — , V x = - - ' x — , (3.19) 

Vpa + Pb WPa + Pb) x WPa + Pb) x 



where x = 1, 2 and the condensates O x and V x have been defined in (3.13). 



J Somc further technical detail is given in Section (4.2) 
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Figure 3: The plots correspond to the two possible ^-condensates in the balanced pa = Pb 
case. The solid black line represent the non interacting setting (7 = 0); red dotted lines 
correspond to increasingly negative (repulsive) values of 7 while green dashed lines refer 
to increasingly positive (attractive) values of 7. The A-condensate is not plotted. We have 
competition/enhancement corresponding respectively to attractive/repulsive interaction. 
In the balanced case the critical temperature is not affected. 

4 Coexistence and mutual competition/enhancement 
of the two condensates 

We want to study how the presence of the A-condensate affects the condensation of ip; in 
particular, we are interested in how 7, the direct coupling between ip and A, influences the 
mutual competition/enhancement of the two order parameters. We plot the ^-condensate 
as a function of the temperature T in the canonical ensemble, namely keeping the charge 
densities pa and pb fixed, in various situations: with and without the A-condensate and 
for different values of 7 (of either signs) . In Figure [2] we show two examples of coexisting 
ip and A condensates. 

In Figure [3] we show the effects of the 7 interaction on the plotted ^-condensates (for 
both the two different quantization schemes) in a balanced case, namely pa = Pb- The 
black solid line represents the 7 = case where the two condensates are completely inde- 
pendent; in this case the presence of a A-condensate has no effect on the ^-condensate. The 
green dashed condensates correspond to increasing attractivd^j] (i.e. 7 > 0) 7-couplings 
while the red dotted condensates correspond to increasingly repulsive interaction. In 
this balanced case the system is symmetric under the exchange of the "electric" and 
"magnetic" sectors and the A-condensates (not plotted) are identical to their ip cousins. 
Attractive values of 7 disfavor the condensates attenuating their amplitudes while neg- 

11 The attractive/repulsive character of the 7 interaction can be checked comparing to the case in which 
the interactions between tp and A are mediated by a would-be scalar (i.e. with even spin) particle. 
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Figure 4: Plots showing the two ^-condensates in the unbalanced case. We observe 
again (as in the balanced cases plotted in Figure |3| competition/enhancement for attrac- 
tive/repulsive interaction but here we have in addition that the critical temperature is 
accordingly affected. The plots are traced for pa ^ Pb with pa and ps kept fixed; the 
temperature is normalized with respect to Tc, i-e. the critical temperature for the 7 = 0, 
non-interacting case. 



ative values of 7 have the opposite effect. Notice that in this balanced case there is no 
effect on the critical temperature at which condensation occurs. 

Turning to an unbalanced case, pa 7^ Pb, we still have the shrinking/amplifying pat- 
tern already occurring in the balanced case but, in addition, also the critical temperature 
at which tp condenses is affected. In particular, attractive interactions lower the critical 
temperature while repulsive interactions raise it. This is in line with the idea of compe- 
tition/enhancement already suggested in the balanced case. Indeed a higher value of the 
critical temperature corresponds to an easier condensation. We show our results for the 
unbalanced case in Figure |4j 

In [9] it was conjectured that two holographic orders are competing or mutually en- 
hancing when their interaction is respectively attractive or repulsive. The present model 
offers then a neat example where this conjecture is confirmed. In Eq. (3.16) and (3.17) 
the mixing terms containing 7 could be thought of as contributing to the effective mass 
of the scalars. As a consequence we have for instance that for ip the presence of a non 
-vanishing A affects its effective mass. We furthermore see that a negative 7 lowers the 
effective mass while a positive 7 has the opposite effect. In relation to the BF bound crite- 
rion for scalar instabilities we know that a lower mass value encourages the condensation 
(see Subsection 4.1). 
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4.1 Comments on stability 

In [7J, studying the unbalanced holographic superconductor with a single order param- 
eter, the criterion for IR stability of a scalar field on the doubly charged ArfS'-Reissner- 



Nordstrorm black hole was given _ The supporting argument is based on the analysis 
of scalar fluctuations on the Reissner-Nordstrorm background in the near-horizon region; 
there the scalar equation turns out to be that of a free scalar on AdS-z geometry with 
appropriate effective mass. Eventually, the AdS% Breitenlohner-Freedman bound gives 
the criterion for IR stability of the scalar. An analogous approach appears to be not 
feasible here. In principle a similar near-horizon analysis could be performed, but here 
we encounter a problem: actually we would like to study one scalar on a hairy black hole 
background (namely on a solution where the other, say A, scalar has already undergone 
condensation); even though we have an analytical expression for the fields of the hairy 



black hole in the IR region 15 , on such background the equation for ip fluctuations does 
not reduce to a simple free scalar on an AdS spacetime. In other words, here we cannot 
just apply a simple BF bound argument to study IR stability. 

Another aspect related to the "stability" of our solutions in a broader sense concerns 
the numerical approach. In fact, the numerical computations performed to study the 
thermodynamics of the holographic system at hand are rather delicate and sensitive. So- 
lutions with coexisting ip and A condensates are sometimes numerically unstable. An 
excessively low numerical accuracy can lead to unphysical "transitions" between a double 
condensate phase to a single condensate phase along the numerical sweep in temperature. 
A numerical instability must however not be confused with a thermodynamical one. In- 
deed, in the holographic framework, the numerical computations consist essentially in the 
numerical solution of a system of coupled differential equations, i.e. the classical equa- 
tions of motion of the dual bulk gravity theory. Said otherwise, we are not performing 
any direct thermodynamical sampling of the partition function of the boundary theory. 

The thermodynamical stability or metastability of the coexisting phase can be checked 
explicitly studying the value of the renormalized Euclidean on-shell action of the corre- 
sponding dual solutions. For the sake of simplicity we have checked successfully the case 
where the scalar condensates have dimension 2 (i.e. we chose the UV boundary conditions 
which set the fastest fall off to zero for both the scalars); in this circumstance the on-shell 



Euclidean action is finite 17 



Throughout this paper we have considered both signs for the quartic coupling 7. One 
could worry about the stability of the scalar potential and particularly about the fact 
that it could be unbounded from below. However, we work in a probe approximation 
where the fields are regarded as small fluctuations; consequently, our effective model can 
be thought of to be corrected by higher order terms in the potential which nevertheless 



See 15 for the analogous study in the balanced case. 
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are negligible as long as the scalar fields are small fluctuation. 



4.2 A technical comment on the code 

In the first model of a probe holographic superconductor [2] , the code employed kept the 
horizon radius (and consequently the temperature) fixed. The sweep in temperature was 
actually performed varying ■0( r i?) an d letting p free; in fact in E) the condensate is plotted 
with respect to the ratio T/p and, taking advantage of a scaling analogous to (3.10), the 
plot is interpreted as a sweep at fixed p and progressively lower temperature. We adopt 
a similar approach here. We fix the horizon radius th to one and require to the shooting 
procedure to keep fixed the ratio between the two charge densities. In this way the plotted 
points corresponding to various values of the ratio T / \JpA + Pb can be interpreted as a 
sweep performed varying the temperature and keeping the charge densities fixed. 



5 Conclusion and future developments 

We introduced a simple holographic model presenting the coexistence of two concomi- 
tant, non-trivial order parameters. The model represents a simple generalization of the 
standard holographic superconductor. Our focus was on the thermodynamical equilib- 
rium properties of the system with particular attention on the mutual interaction of the 
two orderings. Contrarily to the general weak-coupling expectation, multiple orderings at 
strong coupling can coexist and even enhance one another. We observed the possibility 
of such an enhancement and, relying on a probe analysis, we found results in accordance 
with the conjecture advanced in [9] claiming that an attractive/repulsive interaction be- 
tween the dual bulk fields leads to competition/enhancement of the corresponding order 
parameters. Our model possesses two charge densities associated to two Abelian symme- 
tries U(1)a, U(1)b] working in the canonical ensemble we studied the coexistence of two 
orderings both in an unbalanced (pa 7^ Pb) an d in a balanced (pa = Pb) setting. 

Ongoing analysis of the backreacted case indicates that in general, since it is an at- 
tractive interaction, gravity tends to render the two condensations competitive. However, 
still on a preliminary level, the situation could be more complicated: gravitational inter- 
actions are indeed not effectively accountable with a simple shift in the effective mass of 
the two scalars as it happened for the direct 7 interaction in the probe analysis. 

One future point of interest consists in the study of transport properties of the un- 
balanced, backreacted holographic superconductor with two order parameters. According 
to [7], in the normal phase, the linear response of the system to external perturbations 
could be accounted for by a single frequency dependent "mobility" function. However, 
in the condensed or doubly condensed phases, the expectation is that possibly the intro- 
duction of a second mobility function could account for the transport properties of the 
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"normal" and "condensed" part of the holographic strongly coupled system. The system 
introduced in the present paper could offer non trivial checks of this expectation. 

Another ongoing interesting perspective (both in relation to condensed matter and 
QCD) is to consider A to be charged under both gauge fields. In the QCD perspective 
this brings our model closer to model a quark-quark condensate and a chiral condensate. 
In the superconductor sense, instead, this situation could represent the coexistence of two 
superconducting condensates, one s-wave and the other effectively mimicking a p-wave or- 
der parameter; said otherwise, we could possibly have two coexisting Cooper-like pairings 
both in singlet and triplet spin state. Still on a rather speculative level, it could be inter- 
esting to further study whether the present model could furnish a simplified environment 
to study a bi-layered holographic superconductor where the two order parameters would 
represents both superconducting order parameters associated to two different layers. 
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